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ABSTRACT 

By using high-resolution ID hydrodynamical simulations, we investigate the 
effects of purely mechanical feedback from super massive black holes (SMBHs) 
in the evolution of elliptical galaxies for a broad range of feedback efficiencies 
and compare the results to four major observational constraints. In particular, 
we focus on 1) the central black hole to stellar mass ratio of the host galaxy, 
2) the lifetime of the luminous quasar phase, 3) the mass of stars formed in the 
host galaxy within the last Gyr, and 4) the X-ray luminosity of the hot diffuse 
gas. As a result, we try to pin down the most successful range of mechanical 
feedback efficiencies. We find that while low feedback efficiencies result in too 
much growth of the SMBH, high efficiencies totally blow out the hot interstellar 
gas, and the models are characterized by very low thermal X-ray luminosity 
well below the observed range. The net lifetime of the quasar phase is strongly 
coupled to the mass ratio between SMBH and its host galaxy, while the X-ray 
luminosity is generally correlated to the recent star formation within the last 
Gyr. When considering the popularly adopted model of the constant feedback 
efficiency, the feedback energy deposited into the ambient medium should be more 
than 0.01% of the SMBH accretion energy to be consistent with the SMBH mass 
to stellar mass ratio in the local universe. Yet, the X-ray luminosity of the hot 
gas favors about 0.005% of the accretion energy as the mechanical AGN feedback 
energy. We conclude that the purely mechanical feedback mode is unlikely to be 
simultaneously compatible with all four observable tests, even allowing a broad 
range of feedback efficiencies, and that including both radiative and mechanical 
feedback together may be a solution to comply the observational constraints. 
In addition to the adopted observational constraints, our simulations also show 
that the ratio of SMBH growth rate over its current mass and the density and 
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temperature distribution of hot gas can be useful observable diagnostics for AGN 
feedback efficiencies. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: formation — 
galaxies: nuclei — methods: numerical 
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Introduction 



The well established empirical correlations between the mass of supermassive black 
holes (SMBHs) and several properties of their host galaxies are providing new insights 
and also new problems for our understanding of massive galaxy evolution. In particular 
the more mas sive SMBHs are hosted in th e more massive bulges, galaxies, and dark 



2000 



matt e r halos ([Kormendv &: Richstone 



Gebhardt et al 



Haring fc Rix 



SMBH (e.g. 



2000 



Ferrarese 



1995 



2002 



Magorrian et a. 



Yu fc Tremaine 



1998 



2002 



Ferrarese fc Merritt 



Marconi fc Hunt 



2003 



20041). and the sha pe of the light profile is related to the mass of the central 



Graham et al. 



20011). Moreover, the co-evolution of SMBHs and host galaxies 
is also supported by the observational evidence of the proportionality between the galaxy 
star formation r ate (SFR) and the mass accretion rate of the SMBHs ove r a broad range 



of red shifts (e.g. 
200% . 



Boyle fc Terlevich 



1998 



Haiman et al. 



2004 



Miller 



2007 



Shankar et al. 



This tight correlation between the stellar mass of the host spheroid and the SMBH 
mass has been tackled by various theoretical explanations that are based on self-regulation 
processes of feedback effects from the accreting SMBHs. When SMBHs are in an active 
phase (i.e. AGN), feedback via either the mass ejection by winds or jets, or alterna tively 



the emitted radiat i on, regulates th e mass accretion rate and the final SMBH mass (IRees 



1984 



Krolik 



1999 



Heckman 



mechanically c 


riven 


Granato et al. 


2004 



King 



2005 



20081). The feedback can be in the forms of radiatiyely o r 



Silk fc Rees 



Murray et a 



energy transportation near the SMBHs (e.g. 



998 



Fabian 



2005 



1999 



<!ing fc Pounds 



Fabian et al 



Begelrnan fc Nath 



such as photoionization and Compton heating ( e.g. 



20071 ). or of a blast wave (e.g. 



Menci et al. 



2003 



200 



Sazonov et al. 



20061 ). of a turbulent 

T 

5|). o f radiative effects 



2005 



Ciotti fc Ostriker 



20081). 



The idea of feedback from SMBHs also has given rise to various hypotheses that 
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could answer important open questions in galaxy and galaxy cluster evolution. For 
example, significant star formation can be prevented by energy deposited from AGNs, 
which suppresses the supply of cold gas, and the further late star formation after the 



bulk of stars form early (e .g. 



2008 



Cattaneo et al. 



2009 



Croton et al. 



2006 



Kormendv e ; 



galaxies to red galaxies (e.g. ILagos et al. 



al 



Cattaneo et al. 



2007 



Khalatyan et al. 



20091 1 and transforms blue star-forming 



20081 ). The AGN feedback may help us explain 



galaxy downsizin g by making the AGN fe edback processes depend on the mass of dark 



matter halo (e.g. 



Scannapieco et al. 



20051 ). Moreover, the feedba ck may alter energy 



and mass distribution o f intergalactic medium in 



Scannapieco & OhI 


2004; 


Scannapieco et al. 


2005 


2007; 


Chatterjee et al. 


2008 


)• 



galaxy clusters ( 


Cavaliere et al. 


2002 




Thacker et al. 


2006; 


Chandran & Rasera 



Independently of the specific model considered, the problem of the self-regulated 
growth of SMBHs and spheroids can be summarized as two questions about effects from 
AGN feedback. First, self- regulation requires that feedback effects have to be well timed 
responding to the growth of stellar mass. Suppose that there is a long time-lag either 
between the start of star formation and the onset of AGN feedback, or between the 
beginning of the feedback effects and heating surrounding gas, then the pace of star 
formation can be faster than the growth of SMBHs. This case would weaken the tight 
correlation between SMBH mass and bulge mass. Therefore, this timing problem is closed 



relevant to how SMBHs are fuelec 



1983 



ShuU 



1983 



Shlosman et al. 



and when they ignite AGN feedback (ICavaliere et al. 



1990 



Quataert fc Naravan 



2000 



Kawaguchil 



2003 



Naravan fc Quataert 



2007 



2005 



Tan fc Blackman 



Miiller Sanchez et al. 



2005 



Jo gee 



2006 



Konig] 



2006 



Davies et al. 



20091 ) . Second, the impact of AGN feedback on host galaxies 



of SMBHs depends on how the gravitational accretion energy is converted to out-f 



mass and radiation from SMBHs, and is transported to the surrounding medium (iKing 



2003 



Proga 



2005 ; 



Murray et al. 



2005 



Proga 



20071 : 



Proga et al 



2008 



Hardee 



owmg 



20081 ). In 
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some circumstances, the energy conversion might even cause the enhancement of star 



form ation instead of suppressing star 



[e.g. 



Silk 



2005 



Ciotti fc Ostriker 



2007 



• qrmation by a 



Pipino et al 



lowin g shock-induced star formation 



20091 ) ■ Very efficient energy output 



and conversion will tend to inhibit the SMBH growth overly, and very inefficient coupling 
of the outflowing energy to the ambient gas would lead to SMBHs more massive even than 
those observed. This energy conversion problem coupled to the timing problem determines 
the efficiency of AGN feedback and the final properties of SMBHs and their host galaxies. 

Therefore, the proper understanding of the timing problem and the e nergy c onver sion 



ElvisI fl2006h 



problem is critical to understand physics of AGN feedback in detail. As 
emphasizes, the structure and physics of AGNs are entangled together, and it is difficult 
to find a direct connection between feedback from AGNs and other properties of galaxies 
without elucidating the complicated connection of feedback effects onto the host galaxies. 
Even though we can begin to ask questions related to AGN feedback on a large scale such as 



galaxy m ergers, galaxy c. 



feedback 



2006 



Cattaneo et al. 



Monaco et al. 



2007 



usters , and cosmological evolution by using simple 



2005- 



Hopkins et al. 



Menci et al. 



2008 



2005: 



Kawata fc Gibson 



Khalatyan et al. 



2005 



descriptions of 



Thacker et al. 



20081 ). we cannot simplify 



the complex physics entangling the SMBHs and their host galaxies without the better 
understanding of involved phys ics, ranging from the small-scale physics of AGNs to the 



larger-scale physics of galaxies (IBlandford 



20011 ) 



In the literature, feedback has been mainly considered in two aspects: radiative 
feedback and mechanical feedback. Radiative feedback is a consequence of radiative heating 



and radiation pres sure by a strong radiation from AGNs, including p 



Compton heating 


(Ciotti & Ostriker 


2001 


Ciotti & Ostriker 


2007 


). And computing t 



Sazonov et al. 



2004 . 



2005 



lotoionization anc 



Fabian et al. 



2006 



since we directly observe the hard radiation from AGNs that causes the ambient heating: 
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there are few uncertainties in computing this type o: 



associated to winds and 


ets from accreting SIv 


IBHs ( 


Tabor & Binnev 


1993 


Binnev & Tabor 


1995; 


Friaca & Terlevich 


1998; 


Veilleux et al. 


2005 




Sprineel et al. 


2005a; 


Konie] 


2006 





feedback. Yet, mechanical feedback is 



and it can be very important although the strength of the effects is highly uncertain. 
Both feedback modes increase the internal energy of the interstellar medium (ISM) and 
may cause local and/or global outflows. But their properties and efficiencies are different 
because radiative feedback (in an optically thin medium) acts almost instantaneously, 
and is transported without mass loading over the entire galaxy. By any accounting the 
radiative energy emitted by the accreting SMBH exceeds the emitted mechanical energy; 
but radiative feedback couples to the ambient medium far less effectively. So there is a 
trade-off in the effectiveness of the two types of feedback that must be investigated in detail. 

In this paper, we extend th e analysis of the purely mechanical feedback models 



described in 



Ciotti et al. 



( l2009al hereafter Paper I), by comparing the results of simulations 
with different feedback efficiencies against four well-constrained observational properties of 
local elliptical galaxies and their SMBHs. In Paper I, new hydrodynamical evolutionary 
models, combining radiative and mechanical feedback, have b e en de veloped from the 



previous purely radiative models ( ICiotti fc Ostriker 



19971 . 



2001 



20071 ). In particular, it has 



been found that purely mechanical feedback models may not be the complete description of 
AGN feedback in Paper I. Though both feedback modes can be important together, firstly 
we ask in more detail whether popularly used purely mechanical feedback models can be 
accepted to explain important properties of SMBHs and their host galaxies in more detail. 
In practice, (1) as most commonly tested in previous research, we estimate the evolution of 
the ratio between the SMBH mass and the host galaxy stellar mass. This ratio has been a 
key diagnosis to test the formation theory of bulge-dominated galaxies since the correlation 
between two masses is found to be tight in observations. (2) the lifetime of the luminous 
quasar phase in the models is also compared to the observationally estimated lifetime (or 



duty cycle 



13 



see 



Martini 



2004 . for a discussion). These two constraints are mainly governed 



by a small-scale physics around SMBHs. (3) another useful probe of the growing stellar 
mass in elliptical galaxies is the amount of the recentl y formed stellar mass, which is we 



constrained by observations in the rest-frame UV (see 



Yi et al. 



2007 



Kaviraj et al. 



2007 



for a discussion). (4) we finally compare the X-ray thermal luminosity o f diffuse hot gas in 



Pope 



2009, for a 



the models with the luminosity observed for local elliptical galaxies (see 
discussion). If the assumed AGN feedback produces effects which are not consistent with 
these observational tests, we may have to reconsider the adopted AGN feedback models and 
their effects in galaxy evolution. 

We investigate how the four observational tests can be used to constrain the AGN 
feedback models, and how the self-regulation process works depending on the adopted 
models. If there is an appropriate efficiency of the mechanical AGN feedback and it is 
the only possible mode of feedback, then the purely mechanical feedback model must 
be consistent with the four observational constraints. In other words, if any values of 
the feedback efficiency cannot produce acceptable simulation results, we may need to 
consider the more general possibility which includes both radiative feedback and mechanical 
feedback modes, as hinted in Paper I. We also suggest various observable diagnostics which 
are related to the physics of AGN feedback, and which can be used to narrow down the 
feedback models and their efficiency, in addition to the four main tests listed above. 

This paper is organized as follows. In §2 we summarize the main components of two 



^The duty 
weights (e.g. 



cycle is commonly defiii ed as the fraction of time on AGN phase with luminosity 



Ciotti fc Ostrikei 



20011 ). But the net lifetime in this paper is simply defined 
by the net duration of AGN with luminosity above some limit. Therefore, the net lifetime 
used in this paper can be slightly different from the luminosity-weighted AGN duty cycle 



multiplied by the Hubble time. 
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purely mechanical feedback models adopted for the simulations. In §3 and §4 we present the 
general properties of simulation results, and comparisons between models and observational 
constraints are given for a broad range of feedback efficiencies. Other possible tests are 
suggested in §5, while the main results are summarized and discussed in §6. 



2. Models 

2.1. Simulation setup 

Here we summarize the main properties of the purely mechanical feedback models, 
which are described in detail in Paper I. In the Appendix, we summarize the important 
input physics such as mass losses from evolving stars, the various time scales of the problem, 
and the treatment of the purely mechanical and purely radiative feedback; in the present 
section only the formulae required for the discussion are given. 

In accordance to several observational constraints, at the beginning of the simula tion, 



the g alaxy model is described by a spherically symmetric Jaffe model for stellar mass (jJaffe 



19831 ). which is imrnersed i n a dark matter halo, so that the total mass density follows a r 



-2 



profile ( ICiotti et al. 



l2009bl ). The effective radius of the stellar mass distribution is about 
6.91 kpc, and the line-of-sight central velocity dispersion is 260 km/s. The total stellar mass 
M* is 2.87 X 10^^ Mq, and inside the half-mass stellar radius the stellar-to-dark matter 
mass ratio is fixed to unity. The obtained galaxy model obeys the observed Faber- Jackson 
relation and the Fundamental Plane of local ellipticals. In all models, the initial galaxy 
model has the central S MBH with the mass Mbh = O.OOIM^,, approximately following the 



Magorrian et al. 



(119981 ) relationship. This assumption is further discussed in §6. 



We simulate the evolution of gas by solving Eulerian equations of hydrodynamics with 
appropriate mass, momentum, and energy sources (see Appendix). All calculations are 



- 10 - 



conducted with a one-dimensional logarithmic radial grid extending from 2.5 pc to 200 
kpc which has 120 bins. The simulation begins at 2 Gyr, i.e. a redshift of 2; ~ 3.2 for the 
LCDM cosmology with = 0.3, Qx = 0.7, and Hq = 70 km/s/Mpc, and continues until 
14 Gyr. We note that the bulk of the mass of elliptical galaxies is already in place at our 
starting epoch. A variety of heating and cooling processes are included as well as spherically 
symmetric radiative transfer in several bands treated in the Eddington approximation. 



We remark that, as in the previous paper ( ICiotti fc Ostrikerll200ll . |2007| ) and Paper 
1, the galaxy model is initially devoid of gas, and has an outflow boundary condition for 
the last radial grid point. Consequently, the ISM is solely provided by the recycled gas 
produced via stellar evolution, and the conflning effect of the intergalactic medium is not 
considered. For these reasons, the ob tained X-ray luminosi ty should be considered as a 



lower limit of the observed values (see 



Pellegrini et al. 



2009h . 



We simulate an isolated elliptical galaxy where the evolution of gas is not affected 
by an y cosmological effects such as galaxy mergers and satellite accretions ( IHopkins et al. 



20081 ) or just cold gas accretion (IKhalatyan et al 



20081 ). The probability of mergers for 
our tested cases since 2 Gyr may be so small for the con sidered mass range that a passi ve 



evolution of the simulated e lliptical galaxy is still valid (iMaller et al. 



2007 



Drory Sz Alvarez 



2006 



Kang et al. 



20081). Moreover, observations of nearby active merger remnants 
show that the local physical process around the central SMBH such as star formation on 
the nuclear disk and supernovae is more imp ortant than what i s simply describe d in galaxy 



merger simulations with AGN feedback (e.g. 



Tecza et al. 



2000; 



Max et al 



2007, for NGC 



6240), supporting the importance of our approach to include AGN feedback physics in 
detail. We also emphasize that the moderate-luminosity AGN activity at low redshift may 



not be causec 
( iGrogin et al. 



by raergers eyen th ough a bright quasar phase may need mergers as triggers 



2005 



Li et al. 



20081). 
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2.2. Physics of purely mechanical feedback models 



As well known, there are four channels to change the total mass of gas in a galaxy 
y. First, winds from stars leaving the main sequence, giant stars, planetary nebulae, SNe, 
etc. contribute a significant amount of mass in gaseous form as the galaxy evolves: over 
the Hubble time, 10% to 30% of the initial stellar mass is ejected by stars. Second, when 
the gas cools (for spherical models usually in a cold shell), star formation can occur, and 
then the newly formed stars will also contribute back gas with their m ass losses. The 



third channel is accretion on to the central SMBH and the nuclear disk (jShlosman et al 



1990 



Nulsen &: Fabian 



the SMBH ( ITan et al. 



2008 



20001) . In fact, cooled gas can form a circumnuclear disk around 



), and then it can be used to form stars in circumnuclear star 



bursts. The fate of the remainin g gas in the accret ion disk is to produce radiation and 



finally being fed into the SMBH (jPavies et al. 



20071 ). In our simulations, we also model the 



secondary mass loss from circumnuclear stars, which also formed from recyc 



mass loss and SNe) in our models, that con tribute to the SMBH accretion ( iShull 



ed gas (i.e . 



1983 



Padovani &: Matteucci 



1993 



Quataert 



2004 ). Finally, the galactic gas can be blown out 



from the galaxy galactic wind/outfiow. 

Of direct interest for mechanical feedback models are the nuclear winds produced by 
the combined effects from the AGN energy output and the Type-II supernovae exploding 
in the star-forming circumnuclear disk. As assumed in Paper I, cold disks will form 



stars, and gas from the disks can be b 



( ICrenshaw et al 



2003 



Krongold et al 



own out by a strong wind from the central SMBH 



20081). In our prescription of the AGN wind loss, the 



mass loss efficiency of the AGN wind r^^^ is defined to be the ratio of the disk outfiow rate 



^Because our simulation does not consider galaxy mergers and (gas) accretion, a possible 
fifth channel of gas infiow is not included in the models as we already mentioned in the 
introduction 
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to the SMBH accretion rate. In Paper I, we introduced two different prescriptions: 

[A] 




= { I (1) 

4 1 + 0.25/' ^ ^' 

where / = L-^Yi/L-Edd is the bolometric luminosity Lbh of the central SMBH in Eddington 

units, and the maximum wind efficiency rj^ is related to the mechanical feedback efficiency 

(see Equation [3]). Therefore, as / increases up to 2 (see Appendix), 77^ approaches rj^ in the 

model of type B. 



In the purely mechanical feedback models, we include only the p hysics of m ass 



Elvis 



2000 



for the 



momentum, and energy transfer from the nuclear wind to the ISM (see 
description of the probable structure of a quasar), while the effect of a jet is not included. 
In particular, the mechanical energy output of the circumnuclear disk is modeled as 

Ldw = CwMbhC^ + eiiC^(l - /rem,h)^^^, (2) 

where Mbh and Mdh*/T*h are the growth rate of the SMBH mass and the SFR of massive 
stars in the disk, respectively (see Appendix). In the second term of Equation |2l en and 
/rom,h are the feedback efficiency by Type II supernovae and the mass fraction of stellar 
remnants for high-mass stars, respectively. Therefore, when the SFR on the circumnuclear 
disk is high or the accretion rate to the SMBH is high, the total energy output from the 
disk wind is high with the contribution of the nuclear wind modulated by e^. 

Our main concern in this paper is Cw which determines how efficiently the growing 
central SMBH can deposit the mechanical energy into its surrounding ISM. Following Paper 
I, we test two models: 

e^, [A] 

I fBl 
4 1 + 0.25/' ^ ^' 

In the both descriptions, rj^ = 1800e^ so that the velocity of the ci rcumnuclear wind is 



about 10 km/s in accordance with observations of broad-line winds (ICrenshaw et al 



20031). 
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In this paper, we maintain the nomenclature of Paper I, and we refer to purely 
mechanical feedback models of type A and B as MA and MB models, respectively. It is 
finally important to recall that, even though radiative feedback effects are not considered in 
the simulations, Lbh has a fundamental importance in the MB models, cLS 6^ IS dependent 
of the luminosity. 

The popular choice of = 0.005 in type A models has been adopted in the majority 
of previous research. For example, 0.5% of the accretion mass is immediately deposited to 



to radiation ( 


Soltan 


1982 


) in Simula 


;ions of 


Hopkins et al. 


2005; 


Johansson et al. 


2009 





Thacker et al. 



(ISpringel et al 



2005a 



( I2OO6I ) also examine a case 



that corresponds to our model A with = 0.005, which easily heats up gas within a 
few Mpc. Other e xamples are found in i nvestigations of the Suii yaev-ZePdovich effect 



from quasars (e.g. 



evol ution (e.g. 



Chatter iee et al 



Scannapiecq et al. 



Robertson et al. 



2008 



Scannapieco et al. 



20081), downsizing of galaxy 



20051 ) . the evolution of the black hole m ass and bulge mas s 



l2006l ). and the formation of elliptical galaxies (e.g. 



Khalatyan et al. 



20081 ). However, ab initio we do not know of either an empirical basis for the adoption of 



this particular coefficient nor calculations from which it could be derived. 

There have be e n som e trials to pin down an acceptable range of e^. For example 



Scannapieco &: OhI (j2004[ ) investigate the cases corresponding to the MA models with 



0.0025 and 0.01, and conclude that the two values fail to explain quasar luminosity 



function and other struct ure formations. But as t hey note, is observationally not well 
constrained by their tests. 



Levine fc GnedinI ( l2005l ) also examine how a constant mechanical 
feedback efficiency affects the distribution of AGN outflows by testing = 0.001, 0.01, 
and 0.005. They flnd = 0.001 as the best value by comparing their models with 



observationally derived fllling fractions of AGN outflows. As these two examples show. 
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the acceptable range of is still controversial for the MA models. However, the different 
purely mechanical models such as our MB models have not been tested extensively yet, 
even though they appear to be more close to the real AGN feedback process in which 
the feedback efficiency is expected to increase as one approaches the Eddington limit 



( iKurosawa fc Proga 



20091). 



For these reasons, we test eight different values of for both model MA and MB 
ranging from 1 x 10~^ to 5 x 10~^ (see Table [1]). Because our goal is to constrain the 
range of acceptable efficiencies and models with the purely mechanical feedback, the tests 
also include unreasonably low and high values of e^. We note that a model without AGN 
feedback effects was shown completely unreasonable in Paper I, producing an extremely 
massive SMBH. 



3. Results of the MA models 

As we explained in the previous section, the family of the MA models, i.e. the 
mechanical feedback models with a fixed feedback efficiency, has been widely adopted in 
various papers as the simplest prescription of AGN feedback. In Paper I, it is shown how 
these models are characterized by a quite rigid behavior with a sharp transition from 
very strong feedback effects to almost no effects, even though the input physics is quite 
elaborate. Here, we explore other aspects of these models, starting with the commonly used 
= 0.005 (Run 6), and then progressively reducing it from Run 6 to Run 1 and increasing 
from Run 6 to Run 8. 
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3.1. Evolution of the MA models 

Figures [1] and [2] present the time evolution of some global quantities in the MA models. 
In all runs, the peak of the SFR (M*) is reached earlier than the peak of the SMBH mass 
accretion rate (Mbh); in some cases the two peaks are coincident within the limit of the 
simulation time resolution. High feedback efficiencies generally delay an initial star-burst 
and a high Mbh, and stop early further star formation and SMBH mass accretion. As found 
in the simulations from the Run 4 to 8 with high values of e^, the effective time-scale of the 
feedback process is so short that any peaked formation of stars and mass accretion to the 
central SMBH does not occur at all. Meanwhile, the low feedback efficiencies make galaxies 
have extended star formation and SMBH mass accretion with multiple peaks. 

Remarkably, the main difference caused by different feedback efficiencies is discovered 
in the rise in stellar mass AM^, and SMBH mass AMbh- If the initial ratio Mbh/M^, needs 
to be maintained at all times, AMbh/AM* is required to follow the same initial ratio 
~ 10~^. Yet, AM* is always higher than AMbh for = ^ x 10~^, while this trend is 
reversed for high feedback efficiencies, as shown in Figure [2l Therefore, our simulations 
imply extreme fine-tuning of the efficiency in order to maintain the SMBH to stellar 
mass ratio at a constant level. 

Different feedback efficiencies also affect the radial structure of gas. As shown in 
Figures |3] and lU models with the high feedback efficiency are more likely to produce a 
high-velocity outflow at a large distance. For example, at r ~ 20 kpc, the outflow velocity 
is about 800 km/s for Run 6 (e^ = 5 x 10^'^). But around the peak time of the SMBH 
mass accretion, the outflow in Run 2 (e^ = 5 x 10^^) has much lower velocity than 500 
k m/s at the same distance. High feedback efficienc i es enh ance the outflow as pointed out 



by 



Scannapieco fc OhI ( 120041 ) and lLevine fc GnedinI (120051) . even though SFR in their cases 



is not as high as in our simulation (iDavid et al 



20061) 
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Heating by stellar processes and the additional AGN heating alters the temperature 



distribution o 
research (e.g. 



the gas within 



Cattaneo et al. 



abor Sz Binney 



galaxies. As already discussed in numerous previo us 



1993 



Binney Sz Tabor 



1995 



Ciotti &: Ostriker 



1997 



2001 



20071 ). energy deposited by mechanical feedback heats up the surrounding 
ISM, resulting in a core structure of hot gas such as seen in Figures [3] and HI Although the 
high feedback efficiency suppresses further star formation and stellar feedback processes, 
even in these cases the energy and mass supply by the strong AGN wind compensates for 
the lack of stellar feedback processes and finally heats up the gas. 

As shown in Figure El the time evolution of the accretion luminosity is strongly coupled 
to the SMBH mass accretion history. However, we note that the luminosity at t ~ 14 
Gyr is similar even for different values of feedback efficiencies, because the difference in 
the SMBH mass accretion rate is small at late times despite the difference of e^. Lbh is 
generally lower than the Eddington luminosity for both high and low feedback efficiencies. 
Bursts of high Lbh above the Eddington luminosity are found only in the models with low 
feedback efficiencies such as Run 2. As we explain in Appendix, our AGN model permits 
a moderate super-Eddington accretion. The AGN feedback with the low efficiency is not 
strong enough to stop the development of cold gas when the precursor of a large amount of 
cold gas reaches the central region. Therefore, a large fraction of the cold gas is accreted to 
the central SMBH at the end. 

We find that the star formation history is mainly parallel to the evolution of both the 
X-ray luminosity from the hot diffuse ISM and the infrared (IR) emission by dust, which 
re-radiates the absorbed stellar radiation, as shown in Figure [51 As the difference in the 
star formation history implies in Figure (H the difference in X-ray and IR luminosity is more 
significant for different feedback efficiencies than the difference in Lbh- In addition, the 
X-ray luminosity of hot ISM intricately changes more than the IR emission, corresponding 
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to the energy input from the AGN feedback and the mass loss by the continuously escaping 
hot gas (see Figure [2]) • 



3.2. Are the MA models acceptable? 



Figure |6] summarizes the results of our tests on the MA models. The mass ratio 
between the central SMBH and the host galaxy at 14 Gyr, i.e. z ~ 0, monotonically 
increases as the feedback efficiency decreases. For > 1 x 10"'*, the growth of stellar and 
SMBH mass is so small that the deviation from the initial ratio is not significant. Following 
the increased SMBH accretion, the net lifetime of bright AGN phase is also long in the 
models with the lo w e^. When w e adopt the B-band magnitude Mb < —23 mag as a limit 



of optical quasars (IMartini 



(Elvis et al 



(IMartini 



20041 ) and use the typical spectral energy distribution of quasars 



1994[ ). the observational constraint on the maximal net lifetime is about 1 Gyr 



20041 ). Even though there is no well defined observational limit on the net lifetime 



when the bolqmetric luminosity Lbh is higher than 10% of the Eddington luminosity ^Edd, 



Hopkins et al. 



(120051 ) claims that it is comparable to the lifetime of the optical limit on 
quasars. We find that the two measurements can be comparable for only ^ 10"^. But if 
our models have to reproduce any quasar phases since 2 Gyr, ^ 10^^ is not acceptable 
because it does not permit any strong SMBH accretion phases as shown in Figure [1] and 
does not experience any luminous phases with Mb < —23 mag or Lbh > O.lLsdd- 

The episodic quasar lifetime can be used to constrain our models in addition to the net 
quasar lifetime. The episodic lifetime is measured for each instance of the quasar phase, i.e. 
the high SMBH accretion phase, while the net lifetime is the sum of the episodic lifetime 
for all instances. Because our models do not provide any information about quasar activity 
before the initial time of simulations, i.e. 2 Gyr, the net quasar lifetime of the models 
can be increased up to 2 Gyr further. Therefore, the episodic quasar lifetime can be a 
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better diagnostic in our simulations. In Figure [TJ we present the change of the episodic 
hfetime for Lbh > -^Edd as an example. This episodic lifetime does not change significantly 
for every episodic activity, having a typical duration 0.4 Myr which is longer than both 
a fixed simulation output time-step size 0.1 Myr and varying computational time-step 
sizes (in average, about 50 year). The episodic lifetime for Mb < —23 mag also has the 
same pattern, but is abo ut 1 Myr which is acceptable compared to a recent measurement 



( iKirkman fc Tytler 



20081 ) 



We also compare the predicted X-ray luminosity of the hot ISM and the mass fraction 
of recently formed stars to the available observational constraints. The constraint on 
the X-ray luminosity i s derived from the typical X-ray luminosity of local ellipticals 



( jo 'Sullivan et al. 



the X-ray luminosity (jCiotti et al. 



2001) after subtracting the contribution of di screte X-ray sources from 



1991 



Kim fc Fabbiano 



20041 ) . As shown in Figure El if 
the feedback efficiency is too low as in the model with = 1 x 10"^, too much mass is 
accreted to the central SMBH, decreasing the X-ray luminosity and suppressing the late 
star formation. Meanwhile, the low SFR in the models with high feedback efficiencies 
causes the lack of hot gas, which sequentially results in low X-ray luminosities. Although 
the temperature of the hot gas is higher in models with the high feedback efficiencies, the 
total amount of gas is too low to be compensated by the high temperature. However, it is 
important to recall that all models explored in thi s paper represent a ga laxy that is initially 



devoid of gas and without external pressure (see 



Pellegrini et al 



20091 ). The measured 



recently formed stellar r nass for the 



galaxies (lYi et al. 



2007 



Donas et al. 



ast Gyr is lower than a fe w percent in local elliptical 



2007 



Kaviraj et al. 



20071 ). which is much higher than 



what we find in our MA models. Hence, the fraction of recently formed stellar mass in 
simulation results is acceptable compared to local ellipticals. 



The duration of SFR M^, > IMq/ut depends on when intensive star formation occurs 
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and how strong the feedback efficiency is. Generally, the late star formation has a longer 
duration. For example, the late star formation in Run 1 continues longer than 100 Myr, as 
shown in Figure [71 The effect of different feedback efficiencies is significant in modulating 
the early star formation. For example, in Run 3 the peaked early star formation lasts for 
about 10 Myr. Yet, the low feedback efficiency in Run 1 allows the accretion to the central 
SMBH to occur more frequently, resulting in the short duration of vigorous star formation. 

In short, it is difficult to find the range of the feedback efficiency in the MA model 
that satisfies the four observational constraints together. Low efficiencies (e^ < 5 x 10~^) 
produce too massive central SMBHs or long net lifetime of quasar activity. High efficiencies 
(e^ > 5 X 10^^) have different problems: too low X-ray luminosity or no quasar activity. 
Importantly, the popularly used = 0.005 does not pass the four tests simultaneously. 



4. Results of the MB models 



We now move to discuss the MB models. In these models, the mechanical output 
from the central SMBH depends on the accretion luminosity, increasing as a function of 
L-Bn/L-Edd as presented in Equation El This description is definitely more close to the real 
processes ar ound the central SMBHs than the fixed mechanical feedback efficiency in the 



MA models ( iKurosawa &: Proga 



20091 ). As we will see, however, the basic consequence of 
varying the feedback efficiency is qualitatively same in MA and MB models, despite their 
differences. 



4.1. Evolution of the MB models 

The peak value of the feedback efficiency determines how frequently the SMBH and 
its host galaxy can achieve high SMBH mass accretion rate and SFR, as we already found 
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in the MA models. Figure |8] presents the SMBH mass accretion rate and SFR in the MB 
models. Either extremely high or low feedback efficiencies does not permit the resurrection 
of both high Mbh and M^,. For example, in Run 1 SFR is higher than the SMBH mass 
accretion rate at any time, while the reversed pattern of growth rates is found in Run 8. 

The peak feedback efficiency also determines the onset of the earliest burst and the last 
burst in Mbh and M*. In Run 1 with = 1 x 10~^, the initial effect of mechanical feedback 
is too weak to suppress the growth of the SMBH mass, even though star formation is always 
calm. As seen from Run 2 to 7 in Figure [HI increasing prevents the early intensive mass 
accretion onto the BH, but the global SFR is not initially affected by the high feedback 
efficiencies. High feedback efficiencies also cause the early cessation of repeating peaked 
high Mbh and M*. For example, the last peak Mbh of Run 7 is found to be about 5 Gyr 
earlier than that of Run 5. 

The change of the total mass in each run is summarized in Figure [91 We do not find a 
strong variability in the mass of the ejected gas and gas inside a galaxy despite the large 
difference in feedback efficiencies, while we found a significant difference among the MA 
models. The main difference is found in the total mass of stars and the central SMBH, 
following the difference in the evolution of Mbh and M* in Figure [H] In all simulations 
except for Run 1, the increase in stellar mass overwhelms the growth of the central SMBH 
before approximately 8 Gyr. As we find in the MA models, the MB model also shows that 
the dynamically evolving model of AGN feedback naturally results in the time-dependent 
mass ratio between the central SMBH and its host galaxy. 

The importance of physics-based feedback models is obviously found by comparing 
Figure [2] and [91 Even for the same radiation conversion efficiency from the accretion 
mass, the MA model is more effective than the MB model in supplying feedback energy 
and suppressing the growth of both stellar and SMBH mass. Although the maximum 
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instantaneous feedback effect is same for tlie same in both models, tlie liigher feedback 
effect in the MA models at low accretion rates leads to quite important differences. This 
difference finally affects how frequently intensive AGN activity and star formation is 
restored. 

The impact from the differences between MA and MB models is particularly apparent 
for very low feedback efficiencies. For Run 1 with = 1 x 10~^, the SMBH accretion rate 
is almost constant at the high value of 10° '^Mo/?/r in the MB model, while the rate strongly 
fluctuates up to about 8 Gyr in the MA model. The luminosity dependence of e„ in the 
MB models allows more rapid growth of the central SMBH than in the MA models. And 
then the very large value of L^dd with the large Mbh consequently reduces following 
Equation [3] in the MB models. Therefore, the evolution quickly turns into a runaway state, 
while the central SMBH grows. Even though either very low or high values of does not 
show bursting activity (see Figure [8]), in Run 1 Mbh is much higher than that in Run 8 
because the runaway process prevents bursts with increasing Mbh in Run 1, but simply 
high feedback efficiency in Run 8 prevents both the bursts of Mbh and the increase in Mbh- 



4.2. Problems of the MB models 



As we find in the tests with the MA models, the family of the MB models also appears 
to fail to pass our tests. Low feedback efficiencies produce too massive central SMBHs and 
too long net lifetime of quasar activity as presented in Figure [TOl while the ratio M^r/M^ 



is much higher than the locally found value of about 10 . Moreover, the net 



quasar phase at Mb < —23 mag is not matched to the limit of 1 Gyr (iMartini 



ifetime of 



200J). As 



found in the MA model, the recently formed stellar mass in the MB models is also much 
lower than the observational maximum limit. However, the thermal X-ray luminosity for 
any value of except for Run 1, does not conform to the observational limits. 
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The episodic lifetime of quasar activity and significant star formation does not have a 
systematic difference in the MB models compared to the MA models. The duration of the 
last burst is considerably longer than those at earlier times as shown in Figure [TT] The first 
star burst is also longer than other peaks of SFR in the models with low feedback efficiencies 
such as Run 7. The duration of the bright quasar phase weakly changes considering the 
limit of Lboi > -^^Edd- But we note that the duration of episodic activity is subject to how 
we define the limit of activity. 

The family of the MB models also fails to conform to four major observational 
constraints considered together. The fraction of recently formed stellar mass does not play 
a key role as an important diagnosis because of its poor sensitivity to variation of the 
feedback efficiency. The net lifetime of quasar phase is strongly coupled to the change of 
the SMBH mass accretion, as we find in the MA models. The X-ray luminosity of hot ISM 
turns out to be a useful diagnosis to test the feedback physics again in the MB models. The 
X-ray luminosity favors < 5 x 10^^. But this range is much lower than > 10^^ which 
produces the right ratio between the SMBH and its host galaxy mass. 



5. Other possible diagnostics 



In addition to four tests adopted in this paper, we find additional possible diagnostics 
which have not been well investigated observationally, but which can be valuable to 
constrain various AGN feedback models. In particular, the feedback efficiencies need to 
be refined more precisely even in the purely mechanical feedback models because of the 
possibility that the feedback efficiency may be dependent of the SMBH mass and other 



properties (e.g. 
quite important 



Merloni fc Heinz 



20081 ) . Therefore, finding an effective tool of diagnosis is 
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The current ratio of the SMBH accretion rate to the SMBH mass can be a constraint 
on the feedback efficiencies for different feedback models. As shown in Figures [H [HI and [T21 
the cumulative effects of the self-regulati on process result in the differ ence in the SMBH 



mass and the current accretion rate (see 



Netzer fc Trakhtenbrot 



20071 . for a discussion). 



This ratio Mbh/Mbr can be converted to the observable quantity L-Bu/LEdd which are 
lower than 10^"^ in local low-luminosity AGNs and 0.1 to 1 in classical luminous local AGNs 

ni — [ 

( iHol l2008f ) . As we find in Figures [71 and [TH this ratio changes quickly within 1 Myr even 
though the duration of L^n/LEdd > 1 is quite insensitive to the feedback models and their 
efficiencies. The main concern in observation is measuring the statistic al distribution o f this 



ratio for broad ranges of AGN activity levels as well as a quite phase (iKim et al. 



2008h . 



The central den sity and temperature of hot ISM is also an interesting quantity of the 



feedback model (e.g. 



Pellegrini et al. 



20091 ). The simulations show that the total amount of 
gas and its temperature change responding to the energy input by different AGN feedback 
models. In addition to this global change, the change of hot ISM in the central region 
around the SMBH shown in Figure [13] can be compared to observations for the same range 
of the central SMBH and host galaxy mass. For example, the observed central density of 
electrons in a few local quiescen t ellipticals is about 0.02 cm~^ which cor responds to the 



mass density p ~ 10 g/cm ( iHumphrey et al 



2006 



Soria et al. 



2006al lbl). which is not 



very different from the values presented in Figure [131 But as the observational values also 
depend on the specific galaxy models adopted, proper comparisons can be obtained by 
focusing on a well observed galaxy for which the dynamical and structural properties are 
well constrained by independent studies, or by studying statistically a large set of observed 
profiles of density and temperature for the X-ray emitting gas. Because we do not include 
recent accretion of external gas or major/minor galaxy mergers, a direct comparison of 
our simulations results to observations need to be limited to local ellipticals not showing 
dynamical and structural properties of recent accretion or mergers. Moreover, the best 



- 24 - 



comparison between our simulations and observations is possible only with the simultaneous 
measurements of the SMBH mass and its accretion rate too. 



Discussion and conclusion 



We have shown that the implementations of two simple classes of purely mechanical 
feedback models are not likely to satisfy the four major observational constraints when 
considered together. Importantly, there does not seem to be any possible range of the 
mechanical feedback efficiencies, including the commonly used = 0.005, which can be 
compatible with the observations. The simulation results also prove that the self-regulation 
process by AGNs requires a careful consideration of both the timing problem and of 
the energy conversion problem. As we will discuss in a subsequent paper, the models 
including both mechanical and radiative feedback may be a right approach to simulate the 
self-regulation process as we observe in local luminous quasars. 

The simulations presented in this paper do not consider departures from spherical 
symmetry, nor chemical evolution of stellar and gas metallicity in addition to other 
shor tcomings such as the effects of ext ernal gas pressure and accretion of stars to SMBH 



(e.g. 



Jogee 



2006 



Pellegrini et al 



20091 ). Because of these missing parts with the present 
approach, it is not possible to exploit other well-derived observational constraints on the 
feedback models. For instance, the commonly observed extended emission line regions 
around quasars are mainly understood as the consequence of geometrically complex 
outflows driven by eith e r radiative or mechanical wind during th e intensive AGN phase 



( iBatcheldor et al. 



2007 



Let awe et al. 



2008 



Fu fc Stockton 



20091 ). Testing simulations 



based on the properties of the outflows will require one to resolve three-dimensional 
spatial structure, as various hydrodynamic or radiative instabilities are not caught in 
one-dimensional simulations. We note, however, that the adopted prescriptions for the 
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mechanical AGN feedback effects are based on sub-grid physics corresponding to mass, 
momentum, and energy transportation in small-scale turbulent motions, even though 
the simulations assume a spherical symmetry. Implementing chemical evolution models 
in our simulations would be a necessary step to use the observed metalli city gradients 



and 


average i 


netal] 


1993 




Gibson 


1996) 



licities of elliptical galaxies as another diagnostic (e.g. 



CaroUo et al. 



19961 ). Even though the effects from AGN feedback on galactic metallicity 



gradients have not been care 



( jSanchez-Blazquez et al. 



ully considered in the current observational investigations 



20071 ). this assumption needs to be verified in the AGN feedback 



models with chemical evoluti o n by comparing simulati o n results with observed grad ients 



(Matteucci & 


Padovani 


1993 




Friaca & Terlevich 


1998 




Aneeletti & Giannone 


Ballero et al. 


2008; 


Kisaka et al. 


2008; 


Spolaor et al. 


2009 


)• 



2003 



In our simulations, we do not consider the problem of the common evolution of the 
SMBH mass and of the stellar mass during the initial phases of galaxy formation. The 
simulations begin with a galaxy already formed and a central SMBH. In particular, here we 
assume th at the initial SMBH ma y have formed with the right mass ratio for the initial 



bulge (e.g. 



Elmegreen et al. 



20081 ). But the establishment of the SMBH-to-stellar mass 



ratio, as well as the Faber- Jackson and the Fundamental Plane in the early evol ution o f 



Ciotti 



ellipti cals, is still poorly understood both observationally and theoretically (see 
2009d . for a review). It will be interesting to test different initial central black hole to stellar 
mass ratios in future simulations. Because stellar mass determines the total amount of 
recycled gas, and AGN feedback models depend on SMBH mass, probably, different initial 
mass ratios and masses might result in complicated evolution which need to be carefully 
tested in simulations. 

Minor mergers and cosmological gas accretion might not affect our main conclusions. 
As we explained in the introduction, major galaxy mergers might not be a significant effect 
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on late galaxy evolution. But late minor mergers and accretion can play an important role 
in fueling gas onto SMBHs. Even in these cases, our main conclusions might not be strongly 
changed except for the test with the net lifetime of the luminous quasar phase. The net 
lifetime of the quasar phase depends on how many times strong gas fueling occurs, while 
other tests such as the simultaneous growth of SMBH and stellar mass and the episodic 
lifetime of the quasar phase is dominated by a feedback physics. Cosmological simulations 
including detailed AGN feedback effects such as our models will be a direct test of effects 
from minor mergers and gas accretion. 

Despite their intrinsic limitations, our simulations bridge the gap between simple 
prescriptions of AGN feedback in cosmological studies or simulations of galaxy mergers, 
and examination of the small-scale physics around the central SMBHs. The problem of 
the self-regulation process itself demands elaborate prescriptions of dynamical evolution 
and energy conversion around the central SMBH and its surrounding ISM. In this paper, 
we found additional evidence that purely mechanical feedback models, including some 
improved versions based on physical arguments, fail to pass basic observational constraints. 
Our next paper will present how combining both radiative and mechanical AGN feedback 
effects can produce the evolutionary models of elliptical galaxies that are more consistent 
with the properties of the central SMBHs and their host galaxies in local universe. 

We are grateful to Michael Strauss, James Gunn, GiUian Knapp, Renyue Gen, and 
Ghristy Tremonti for useful discussions and careful reading. We thank the anonymous 
referee for considered comments which improved this manuscript. M.-S. is supported by 
the Gharlotte Elizabeth Procter Fellowship of Princeton University. Computations were 
performed on the computational facilities of PIGSciE (Princeton Institute for Gomputational 
Science and Engineering). 
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A. Input physics 



In this Appendix we summarize the implementation of the physi cs involved in the 
simulations. For a more extensive discussion the reader is referred to 



Ciotti &: Ostriker 



(120071 ) and to Paper I. 



A.l. The hydrodynamical equations 

The evolution of the galactic gas flow is obtained integrating the time-dependent 
Eulerian equations of hydrodynamics: 

l^ + V- M = ap, + pii-p+ (Al) 
— + V • {mv) = -(7 - l)VE - Vprad + 9P- m+ (A2) 

— + V-{Ev) = -{-f-l)EV-v + H-C+ (A3) 
at 

jap* + Pu)iv^ + Sal) , p , p p+ 

p, m, and E are the gas mass, momentum and internal energy per unit volume, respectively, 
and V is the gas velocity. The ratio of the specific heats is 7 = 5/3, and g{r) is the 
gravitational field of the galaxy (stars and dark matter), plus the contribution of the central 
SMBH. The gravitational field is updated at each time step by considering the SMBH mass 
growth; for simplicity, we do not take into account neither the ISM contribution, nor the 
mass redistribution due to the stellar mass losses and star formation. The total radiative 
pressure gradient is Vp^ad = (Vprad)cs + (Vj9rad)dust + (Vprad)photo, while if - C is the 
radiative heating and cooling term. 

The energy source term is obtained under the assumption that the streaming velocity 
of the source distribution is zero, neglecting the small contributions of the internal energy 
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of the injected gas, and of the kinetic energy of stellar wind when compared to the 
local stellar velocity dispersion contribution (for the derivation and detailed discussion 
of the hydrodynamic al equations with moving isotropic or anisotropic source terms, see 



D'Ercole et al. 



( l2000l )). The source terms ap* and Ei of the initial, passively evolving stellar 
population, and the source terms due to Type II Supernovae, pn and En, are described in 
the following. 

The first active grid point Ri is placed within the Compton radius 

2GMBHpmp Mbh lO^K 

so that at Ri we can impose the physical condition of a vanishing thermodynamical pressure 
gradient, leading to gas free-fall on the circumnuclear disk when the radiation pressure is 
negligible; in this paper we adopt Tx = 2.5 x 10^ K. The appropriate values for radiation 
pressure at Ri are obtained from the disk treatment. 

The simulations are realized with a spatially second-order Eulerian scheme which 
adopts two staggered grids, each of them consisting of 120 logarithmically spaced grid 
points. The equations are integrated with a time-splitting scheme, while the heating and 
cooling terms in the energy equation are integrated by using a predictor-corrector scheme, 
so that the integration is second order in time. At each simulation time, the time-step is 
determined as a fraction of the minimum among the Courant condition over the grid, and 
of the others characteristic times associated with the described physical processes: during 
the accretion phases (and subsequents bursts of radiation), it is not infrequent to have 
time-steps of the order of 1 yr or less. However, it is important to note the accretion events 
are characterized by the intrinsic time-scale related to equation (1A4I) by 

where cx is the isothermal sound velocity associated with the Compton temperature. 



= 1^ ~ 1.22 10 V'/'^^ ^] yr, (A5) 
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A.2. Stellar peissive evolution: SNIa rate and stellar mass losses 

The stellar mass loss rate and the SNIa rate associated with the initial stellar 
distribution are the main ingredients driving evolution of the models. In the code the 
stellar mass losses - the source of fuel for the activity of the SMBH - follow the detailed 
prescriptions of the stellar evolution theory. Over the whole galaxy 

= IMF(Mto)|Mto|AM, (A6) 

where the initial mass function IMF is a Salpcter law (normalized as described in CDPR), 
and the turn-off mass (in Mq) of stars at time t (in Gyrs) is 

logMTo = 0.0558(logt)2 - 1.338 log t + 7.764. (A7) 

Finally 

r Mto - Mfin(MTo) = 0.945Mto - 0.503, (Mto < QMq), 
AM = { (AS) 
1^ AM = Mto - 1.4M0, (Mto > 9M0). 

The time evolution of the SNIa rate is parametrized as 

-12r2n Lb f t \ ' , 



R^^it) = 0.32 X 10-"/.^^s«^ ) yr-'. (A9) 

where h = i^o/lOO km s~^ Mpc~^, and the coefficient i^sn allow for different choices in the 
present-day SNIa. Assuming for each supernova event an energy release of £^sn = 10^^ erg, 
a fraction 77SN of which is thermalized in the surrounding ISM, the energy input per unit 
time over all the galaxy body is given by 

LsN(i) = 1.015 X lO^i/i^sN^N^ (TTVTr-) ^^S^"'; (^^0) 

Lb0 V13.7Gyry 

in this paper we restrict to the case -dsN — 1 and h = 0.75. Here we restrict to the currently 
favoured s = 1.1 value. 

Besides energy, supernovae provide also mass. We assume that each SNIa ejects IAMq 
of material in the ISM, so that the total rate of mass return from the aging initial stellar 
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population at each place in the galaxy is 

^ = (a* + q;sn)p*, (All) 

where q;sn(^) = 1-4Mq i?sN(^)/-^* and a^{t) — M^{t)/M^ are the specific mass return rates. 
With these definitions, the SNIa kinetic energy injection per unit volume in the ISM can be 
written as 

El = r^sN^^SN-^P* = VsnEsn ^ ■ (^12) 

A. 3. Star formation, SNII heating and stEirburst properties 

Star formation cannot be avoided when cool gas accumulates in the central regions of 
elliptical galaxies. In particular, we compute the star formation rate at each radius r from 
the equation 

Pt = 7-form = max(Tcool,Tdyn), (Al3) 

^form 

where p is the local gas density, reform = 0.03 — 0.4, and the associated characteristic times 
are 



_ ^ • / N _ / 3 _ 27rr , . 

Tcool = 77, Tdyn = nim(rjeans,TrotJ, Tjeans = A/ on ^ ^Vot = — TT"- (AMj 

G V oZnup Vc{r) 

E and C are the gas internal energy and the effective cooling per unit volume, while Vc{r) 
is the galaxy rotational velocity at radius r. In the code the stars are maintained in the 
place where they form, and in each shell the associated sinks of momentum and internal 
energy per unit volume are given by the negative of 



^form 

where m is the specific momentum of the ISM. 



^form ^form 
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For a total mass AM,,, of newly formed stars in a given time-step and at a given place, 
we assume a Salpeter IMF 



dM ^^-^'lltj -M^'^WJ ' ^->^<M>M,., = 0.1M^J, (A16) 
SO that the associated total number of Type II Supernovae is 

iV„ = f ^dM = f 1 - f . 7 X 10-^. (A17) 

where the numerical value holds for x = 1.35. As for SNIa, we assume that each SNII event 
releases i?sN = 10^^ erg of kinetic energy, and the resulting mean efficiency is 

in this paper we assume ?7sn = 0.85. The characteristic time for SNII explosion is fixed to 
Tii = 2 X 10^ yr, and from equations (1A13P and (lAlSp their luminosity (per unit volume) at 
each radius from the galaxy center is 

En(t) = ^ r pt{t')e-^'~''^'^^^dt'. (A19) 

ni Jo 

We assume that each explosion leaves a neutron stars of 1.4M0. As a consequence, the 
total mass ejected by the SNII explosions per unit mass is 



1-4-^ ^ 0.2, (A20) 



AM, V Mil / AM, 
and the mass return rate per unit volume of the young evolving stellar population is given 
by 



puit) ^ — I p+(t')e-(*-*')/^Mt'. (A21) 



ni Jo 

Finally, in the code we also compute the fiducial optical and UV luminosity per unit 



volume of the new stars as 

- / 

'opt Jo 



KAt) = ^ f pt{t')e-^'-''^'^°^'dt\ (A22) 

Topt Jo 
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and 

^uv(t) f pt{t')e-^'-n/^-dt', (A23) 

TUV Jo 

respectively, where eopt = 1-24 x 10~^, euv = 8-65 x 10~^, Topt = 1.54 x 10^ yr, and 
Tuv — 2.57 X 10^ yr are the efficiency and characteristic time of optical and UV emission, 
respectively. 



A. 4. Radiative heating and cooling 

Compton heating and cooling, bremsstrahlung losses, line and recombination 
continuum heating and cooling, are taken into account. 

A good approximation to the net gas energy change rate E, valid for T>10'^ K (all 
quantities are expressed in cgs system) is given by 

E^n\Si + S2 + Ss) = H -C, (A24) 

where n is the Hydrogen density (in number), and positive and negative terms are grouped 
together in the heating (H) and cooling (C) functions. The bremsstrahlung losses are given 

by 

Si = -3.8 X IQ-^'^Vf, (A25) 
the Compton heating and cooling is given by 

,52 = 4.1 X lO-^^(Tx-T)^, (A26) 

where Tx is the Compton temperature, and finally the sum of photoionization heating, line 
and recombination continuum cooling is 
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where 



1^ 



80 



17 



g25(logT-4.35)2 g5.5(logT-5.2)2 g3.6(logT-6.5)2 

b=1.7x 10^ T-°-^ 

1.1 4x10^^ 

c = 1.1 — 



,T/1.8105 



2^4 



(A28) 
(A29) 
(A30) 



and 



1 



1.5T-0-5 + 1.5 X 1012 T-2-5 



+ 



4 X 10^° 

y2 



1 + 



80 



,(T-104)/1.5103 



(A31) 



Equations ( IA26P - flA27p depend on the ionization parameter 

eff 



-^BH.photol'^' 



n(r)r2 ' 

where -^^BH,photo(^) the effective accretion luminosity at r, which is evaluated by 
numerically solving in each shell the balance equation 



(A32) 



^-^BH,photo(^) 

dr 



(A33) 



with central boundary condition i^gijpb 



oto V 



0) = LBuit) given by equation flA46p . The 



photoionization+Compton opacity associated with radiation absorption is then obtained 



1 



K 



photo 



'^-^BH,photo('^) 



Anr'^Hir) 



pD 



cS ( 
BH, photo V 



dr 



.photo \ 



Finally, the bolometric ISM luminosity is obtained from equation f lA24p as 



LJr) = 47r / Cr^dr. 



(A34) 



(A35) 



A. 5. Radiation pressure 

Radiation pressure due to electron scattering (where neither the photon numbers, nor 
their energy change) is computed as 

(Vp,ad)cs = : T—, , (A36) 
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where Kes = 0.35 in c.g.s. units, and from equations f lA22|) - flA23P 



L^yir) = 47r / E^yr^dr, Lopt(r) = An Eoptr^dr. (A37) 
Jo Jo 

Note that all the luminosities used in equation (1A36P are unabsorbed. 
The radiation pressure contribution due to dust opacity is given by 

ivpi.adjdust = — 2 — 2 — T' ^^"^^^ 

c 47rr^ c 47rr^ c 47rr^ 



where 



LiR(r) ^ L^'j^,uv(^) + CH,opt(0 + ^uv(0 + Llllir), (A39) 



is the infrared luminosity due to recycling of photons absorbed from the ISM, and we adopt 
as estimates for (cgs) opacity in three bands 

'^opt = Yip^Tio^' " '^^^ " 150 ■ ^ ^ 



At variance with electron scattering the effective luminosities appearing in equations ( ]A38I) 



(1A39P take into account absorption, and are obtained by numerically solving the two 



lowest spherically symmetric moment equations of radiative transfer in the Eddington 
approximation (e.g., Chandrasekhar 1960): 

^ = Anr^E^y - K^ypJtj%), = Anr\E,,, - KoptpJoft)- (A41) 

'^'^uv _ SkuvP-^^uv '^'^opt _ S^optP-^^opt (A42) 



dr 47rr^ ' dr Anr"^ 

— T , V , 

^optV 



The central boundary conditions for stellar luminosities are L^(0) = Ld.uv, -^opt(O) = -^d,opt, 



Juv(O) = Ld,uv/16vr^-Ri and Jopt(O) = -^d,opt/16vr^i?p The effective accretion luminosities 
-^BHUV -^BHopt computed with two equations similar to ( 1A41I) . where the 
distributed source term is missing, J = L^^/Anr'^, and in the UV and optical bands 
^BH,uv(0) = 0.2LBH(t) and Lfj„pt(0) = O.ILbhI^), respectively. 

The last contribution to radiation pressure comes from photoionization opacity, 

/V7 \ P^^photo -^BH,photo('^) /A/IQ\ 
(Vprad)photo = Z 5 • (A43) 
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A. 6. The circumnuclear disk and the SMBH accretion luminosity 

The circumnuclear disk, which is the repository of the gas inflowing at a rate Mf^ 
from the flrst active mesh point Ri of the hydrodynamical grid, and which feeds the central 
SMBH at a rate Mbh? contains at any time the mass gas M^g and a total stellar mass 
Md* = Mdi* + Mdh*, which is divided among low and high mass stars (with the division 
mass at 8Mq). The disk also contains a mass Mrem of remnants from the earlier generations 
of evolved stars. 

In the adopted scheme the accretion rate on the central SMBH is given by 

Mbh = (A44) 



where 

Mm = — Vd = — V — , ^Edd = 7 (A45) 

Td 2Mem eoc^ 

are the flducial depletion rate of gas from the circumnuclear disk, its normalized value, and 
the Eddington mass accretion rate, respectively. The reference radiative efficiency eo is 
defined in eq. flA47l) . Equations flA44l) - flA45P are designed to guarantee that when rj^ ^ I 
the gas is accreted onto the central SMBH at the rate Mm, while Mbh = 2MEdd for r^d ^ 1 
(i.e., we allow for possible moderate super-Eddington accretion; note however that outside 
the first grid point -Ri the fiow accretion rate is limited in a self-consistent way by feedback 
effects). From eq. flA44p we calculate the instantaneous bolometric accretion luminosity as 

i^BH = cem Mbh c^ (A46) 

where 

Am . _ Mbh 

and A is a free parameter so that eEM ~ eoArn for m ^ A~^. In our simulations we fix 
A = 100, and we introduce the normalized accretion luminosity 

/ - ^ = (A48) 
LEdd l + Am' 
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where the last expression derives from the ADAF phenomenological description. 

There are a few lag times in our problem which are expressed as follows. The first is 
the instantaneous disk lag time, appearing in eq. (1A45P : 



where a ~ 10^^ — 10^^ is the disk viscosity coefficient, and -Rd and Mbh are the 
instantaneous values of the fiducial radius of the circumnuclear disk and the mass of the 
central SMBH. We use the scaling predicted by thin-disk theory 

/ M \^^^ 

R,{t) = f,Ri X , (A50) 

where Mbho is the central SMBH mass at the beginning of the simulation. We assume 
fd = 0.4, so that Rd{0) ~ 2 pc for an initial SMBH mass of ~ IO^Mq. 

The second characteristic time is the instantaneous infall lag time from Ri to the disk: 



n = —, vs = \ — , A51 

so that the effective rate at which gas accretes on the disk is obtained by solving the 
differential equation 

mf _ M, - _ ^^^^^ 



dt T\ „ 



where Mi is the instantaneous rate at which gas flows through the flrst active grid point] 
It follows that when Mi provided by hydrodynamics drops to zero the circumnuclear disk 
experiences a fueling declining exponentially with time. 

The disk total gas mass Mdg is not only the source of SMBH accretion, but also of star 
formation in the disk: we assume that a fraction of Mdg is converted into stars at a rate 



^Mi is taken positive in case of accretion and zero in case of outflow at Ri. 
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rj*Mfid (where r^* ~ lOM^g/MBn), and that another fraction of Mdg is lost as a disk wind 
and as a jet at instantaneous rates given by t/^Mbh and r/jMeH, so that the equation for the 
gas mass in the disk is 

^ = Mf -(! + //„ + r/j)MBH - r^^Mfid. (A53) 

The stars formed in the disk are described separately as a function of their mass, i.e., 
high-mass stars (M > Mn = 8Mq) produce a total disk mass Mdh*, and low-mass stars 
(Minf < M < Mil) contribute to a disk mass Mdi* according to the equations 

— — = (1 - fh)V*MM ; —rr- = fhV*MM ■ (A54) 

at T^i at r*h 

For the characteristic evolutionary times we adopt r^,i = Topt and r^,h = th, while we 
assume /h = 0.5, corresponding to a top-heavy Salpeter-like initial mass function of slope 
X ~ 1.16 and minimum mass M^f = O.IM0. The associated optical {Ld,opt) and UV 
(LdjTv) luminosit i es of the stellar disk are calculated following the scheme described in 



Ciotti fc Ostrikeii (120071 ). Finally stellar remnants mass in the disk evolves as 



dMrera_ Mil* ^ . Mih* 

J, — /rcm,l + /rcm,h ; (^AJdJDJ 

at r*i r^,h 

where frem,l = 0.2, /i.em,h = 0.09. 

The equation for the mass loss associated with the disk wind is 

— — = ?7„Mbh + (1 - /rem,l) h (1 - /rem,h) : (A56) 

at r*i r^,h 

the first term is a mass loss driven as a wind by the central SMBH, and the second and 
third are from high mass and low mass stars in the central disk. 

We explore two different classes of models, that we call Type A and Type B, with 

[A] 




[Rl ^^^^^ 

4 1 + 0.25/' 
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We also consider another mass component ejected by disk, i.e. a nuclear jet with 
instantaneous mass flow 

SO that the mass ejected by the jet is always negligible with respect to the wind mass loss 
in Type A models, while it is slightly dominant over the wind in Type B models at low 
luminosity ratios. In the code, all the equations presented in this Section are integrated 
numerically with a first order finite difference scheme. 



A. 7. The mechanical feedback treatment 

We now discuss how the kinetic energy, momentum and mass of the BLR wind are 
transferred to the ISM. The fiducial instantaneous mechanical luminosity of the disk wind 
is given by 

Ldw = e„MBHC^ + eiic^(l - /rom,h)^^^^, (A59) 

where is the mechanical efficiency of the wind, and the second term describes the 
energetic associated with the SNII explosions of the high-mass stars in the circumnuclear 
disk. In analogy with eq. (IA57p . we assume 



e^, [A] 
= <; 36^^ / (A60) 

4 1 + 0.25/' ^ ^' 

In Type A models we explore the range 3 10~^ < < 5 10~^. In Type B models, where 
the wind efficiency is a function of the normalized accretion luminosity, is the maximum 
possible value (reached for 1 = 2). In both cases the instantaneous disk wind velocity is 
given by 

Mdw ^ n 
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where the last expression neglects the mass return contribution of massive stars in the 
circumnuclear disk. In Type A models is in the range 2 x 10^ — 2 x 10^ km s~^ (as a 
function of the specific assumed value for e^), in agreement with observations of BLRs. For 
the same reasons, in Type B models we require — 10^ km s~^, so that 77^ and are 
linked by the relation 

= 1800e^. (A62) 

In analogy with the wind component, the instantaneous jet mechanical luminosity is 
written as 

,v 2 0.0125 
L, = s^Mbkc , ej = — , (A63) 

and the jet velocity is given by 



^j^W— = J— (A64) 

which, for our chosen parameterization gives high but subrelativistic jet velocity of 
v^/c ~ 10~^-^^ for Z^O.l. Finally, the wind and jet momentum are defined as 

mj = Mj'Uj; mw = M^^v^. (A65) 

We now illustrate how we distribute the mechanical feedback over the galaxy ISM. First we 
introduce the instantaneous wind and jet lag times 

Twj = — (A66) 

from the center to the first active grid point R\ (where the subscript indicates the specific 
component - disk wind or nuclear jet - considered), and at each time step we compute the 
time-lagged values for mass, momentum, and kinetic energy at Ri by solving the differential 
equation 

dXi X — Xi , . _. 

dt Twj 
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where Xi is the generic lagged variable associated with the instantaneous unlagged value X. 
Outside Ri we then distribute mass, momentum and kinetic energy over the hydrodynamical 
grid (outside by integrating numerically the phenomenological differential equation 



dlnr Pwjir) v„j dt 



(A68) 



where Y„j is the mass, momentum and energy of the disk wind/jet component at distance r 
from the center, Pwj{r) is the local wind/jet pressure, and for each quantity Y{Ri) = Xi. 
In this paper we restrict to simulations where the time derivative is neglected. In practice, 
we first integrate eq. f lA68P for the wind/jet pressure, i.e.. 



where Y^j is the effective wind/jet momentum crossing the shell of radius r, so that 
eq. f lA68p is a non-linear differential equation for Y^j. Once the equation is integrated, the 
radial behavior of Pwj and the r.h.s. of eq. ( ]A68p are known over the whole grid, and the 



equation can be integrated for mass and energy. 

The solid angle in the denominator of eq. (IA69I1 is the opening angle of the wind and 
of the jet, and the factor of 2 accounts for the biconical nature of the flow. While for the 
jet we assumed in all the simulations the fiducial value Af2j = 2.5 10^^, for the wind case 
we adopt 

(n [A] 

= { , (A70) 

[7rmin(V/2 + a2,l), [B], 

where case B is designed to mimic the behavior found in radiation driven winds: higher 
luminosity corresponds to a larger opening angles. The constant inside the square root is 
fixed to a = Afij/vr, so that for small values of accretion luminosity the wind opening angle 
coincides with the jet opening angle. Finally, note that the almost linear dependence of 
on I for / > a assumes that the linear opening angle depends on \/l for this regime. 



-41 - 



To implement numerically the mechanical feedback terms, we finally compute the 
nuclear wind mass, momentum and kinetic energy per unit volume deposited in each shell 
as 

Source^, = ^3-—^^ (ATI) 

and we add them (only for the wind component) to the r.h.s. of eqs. flAip -(A3). 
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Table 1: Efficiency parameter of the mechanical feedback in the both MA and MB computed 



models. 



Run 






Run 






1 


1 X 10" 


-5 


5 


1 X 10" 


-3 


2 


5 X 10" 


-5 


6 


5 X 10" 


-3 


3 


1 X 10" 


-4 


7 


1 X 10" 


-2 


4 


5 X 10- 


-4 


8 


5 X 10- 


-2 
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t (Gyr) 



Fig. 1. — SMBH mass accretion rate Mbh {thick line) and the SFR {thin line) in the MA 
models. Both rates show a strong time-dependence. The models with the lower feedback 
efficiency generally produces higher Mbh and more frequently. The Eddington accretion 
rate {short-long dashed line) is usually higher than Mbh except for the peak activities. 
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Fig. 2. — Time evolution of the mass budget in the MA models: the change of the SMBH 
mass AMbh {thick solid line), the change of the stellar mass AM* {thin solid line), the 
total mass of gas within 10 Re in the galaxy Mg {thick short-long dashed line), and the 
time- integrated mass of blown-out gas AM^ {thin short-long dashed line). Following the 
difference in Mbh and M*, models with different efficiencies result in different evolution of 
the mass ratio between the central SMBH and stars. Mg decreases because the total SMBH 
accretion rate, star formation rate, and galactic wind rate dominate at late times over the 
stellar mass loss rate. 
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Fig. 3. — Radial structure of the gas internal energy density, velocity, mass density, and 
temperature in the MA model with = 5 x 10~^ (Run 2), at t = 3.5 Gyr {thick line) and 
5.5 Gyr [thin line), corresponding to inter-burst phases (see Figured]). The distributions of 
internal energy density and mass density do not change significantly between the two epochs. 
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Fig. 4. — Radial structure of gas properties in the MA model with = 5 x 10"'^ (Run 6) at 
the same epochs as Figure [31 The comparison with Run 2 (Figure [3]) shows that the velocity 
field and temperature distribution strongly depend on the feedback efficiencies. The high 
feedback efficiency induces the high-velocity outflow, the increase in gas temperature which 
is characterized by a flat high-temperature region of about 100 pc, and the corresponding 
low value of the central density. This trend is consistent with the overall suppression of both 
mass accretion to the SMBH and star formation, as seen in Figure [21 



- 55 - 





Fig. 5. — Time evolution of the dust-extincted optical accretion luminosity Lopt (left panel) 
and the X-ray luminosity from hot ISM within 10 and infrared re-emission of stellar 
radiation (right panel) for Run 2 and 6. Lopt {thick solid line) is basically coupled to the 
mass accretion history shown in Figure [H but it is usually much lower than the 10% of the 
Eddington luminosity [thin solid line). Only models with low feedback efficiencies achieve 
the optical luminosity which is higher than the 10% of the Eddington luminosity for short 
periods. The overall change of the X-ray luminosity from hot gas {thick line) and the infrared 
re-radiation of stellar light by dust {thin line) follows the evolution of SFR rather than the 
history of the SMBH mass accretion. But right after the peaks of the SMBH accretion rate 
in Run 2, the X-ray luminosity {thick line) declines more quickly than the infrared emission. 
This difference implies that the X-ray luminosity may be a more sensitive probe to measure 
the effects of AGN feedback on the ISM than IR emission. 
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Fig. 6. — Results of the four proposed tests for the MA models at the epoch of 14 Gyr. 
From top to bottom, four tests arc the SMBH mass to stellar mass ratio, the net quasar 
lifetime, the X-ray luminosity of the diffuse hot gas inner 10 R^, and the fraction of stellar 
mass formed within the last 1 Gyr fusp- The ratio Mbb./M^ in simulations can be consistent 
with observations for e^^ > 1 x lO""^. Yet, Lx is inconsistent for e^^ > 1 x 10~^, and the 
corresponding models do not experience quasar-like phases. The arrow bars represent the 
range of that is acceptable for each test. For the net quasar lifetime, » 10~^, which 
is represented by the dotted arrow bar, is acceptable only if our models do not need to 
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Fig. 7. — Episodic lifetime of intensive star formation and quasar activity in the MA models. 
For Run 1, 2, and 3, we measure the episodic hfetime of SFR Mq > 1 Mq/ut [thin line) and 
when the bolometric luminosity from the central SMBH Lbh is higher than the Eddington 
luminosity I/Edd {thick line). The typical episodic lifetime of Lbh > ^Edd is about 0.4 Myr 
for all three simulations at any time. The high SFR generally shows a long duration at late 
time. 
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Fig. 8. — SMBH mass accretion rate Mbh {thick line) and SFR (thin line) in the MB 
models. The Eddington accretion rate (short-long dashed line) is always much higher than 
Mbh except for the peak activities, as found in the MA models. In the models with either 
extremely high or low feedback efficiencies, self-regulated outbursts are not found, showing 
almost constant or decreasing Mbh and M*. Compared with the MA models, the MB models 
substantially increase the number of bursts because is low at low accretion luminosity, 
and this favors a quick accumulation of recycled gas in the galaxy. 
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Fig. 9. — Time evolution of the mass budget in the MB models: the change of the SMBH 
mass AMbh {thick solid line), the change of the stellar mass AM^, {thin solid line), the 
total mass of gas within lORe in the galaxy Mg {thick short-long dashed line), and the time- 
integrated mass of blown-out gas AM^ {thin short-long dashed line). Only AMbh depends 
strongly on the feedback efficiency. Because SFR is generally higher than Mbh as shown 
in Figure [H the increase in stellar mass outpaces that of the SMBH mass during the early 
evolution. 
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Fig. 10. — Results of the four tests for the MB models at the epoch of 14 Gyr. Prom top 
to bottom, simulation results are compared to four observational constraints: the SMBH 
mass to stellar mass ratio, the net quasar lifetime, the X-ray luminosity of diffuse hot gas, 
and the fraction of stellar mass formed within the last 1 Gyr. The low in the MB 
models enhances the growth of the central SMBH, resulting in the high values of Mbh/^^* 
and net quasar lifetime. The arrow bars represent the range of that is consistent with 
observations. This comparison shows that none of the tested feedback efficiencies can pass 
all tests simultaneously. 
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Fig. 11. — Episodic lifetime of intensive star formation and quasar activity in the MB 
models. For the same hmits of SFR and the bolometric luminosity from the central SMBH 
used in Figure U\ the MB model also shows the same pattern as the MA model shows. Late 
bursts of star formation {thin line) are maintained longer than early bursts. The duration 
of Lbh > -^Edd {thick line) is between 0.1 Myr to 1 Myr generally without a significant 
dependence on e^. 
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Fig. 12. — Distribution of the SMBH accretion rate versus its mass at 14 Gyr for all 
models in Table [H The numbers near the points are the names of the simulation runs. The 
difference between MA and MB models results is apparent in the different ranges of Mbh 
for the same initial SMBH mass. In both models, high accretion rates correspond to low 
feedback efficiencies. But we note that the central SMBHs in all models are not in an active 
phase at 14 Gyr as shown in Figures [1] and [HI 
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Fig. 13. — Radial temperature and density distribution of hot gas in Run 1, 4, and 7 at 14 
Gyr when the SMBH accretion rate is much lower than the Eddington accretion rate. In 
both MA and MB models, the high feedback efficiencies result in the low central density 
{thin line). However, the temperature profile {thick line) is dependent of which feedback 
prescription is used. In particular, the temperature core in the MA models is produced by 
the constant feedback efficiency, while in MB models the low accretion luminosity at 14 Gyr 
produces a very weak feedback which makes the temperature profile be similar to a standard 
low-luminosity hot accretion profile. 



